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ABSTRACT 

Apart from the 11-year solar cycle, another periodicity around 155-160 days 
was discovered during solar cycle 21 in high energy solar flares, and its presence in 
sunspot areas and strong magnetic flux has been also reported. This periodicity 
has an elusive and enigmatic character, since it usually appears only near the 
maxima of solar cycles, and seems to be related with a periodic emergence of 
strong magnetic flux at the solar surface. Therefore, it is probably connected with 
the tachocline, a thin layer located near the base of the solar convection zone, 
where strong dynamo magnetic field is stored. We study the dynamics of Rossby 
waves in the tachocline in the presence of a toroidal magnetic field and latitudinal 
differential rotation. Our analysis shows that the magnetic Rossby waves are 
generally unstable and that the growth rates are sensitive to the magnetic field 
strength and to the latitudinal differential rotation parameters. Variation of 
the differential rotation and the magnetic field strength throughout the solar 
cycle enhance the growth rate of a particular harmonic in the upper part of the 
tachocline around the maximum of the solar cycle. This harmonic is symmetric 
with respect to the equator and has a period of 155-160 days. A rapid increase 
of the wave amplitude could give place to a magnetic flux emergence leading to 
observed periodicities in solar activity indicators related with magnetic flux. 

Subject headings: Sun: oscillations — Physical Data and Processes: magnetic 
fields — MHD — waves 
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Introduction 



D uring solar cycle 2 1, a short peri odicity between 152-158 days was discovered in 



flares (jRieger et al 



Kile fc Clivei 



electrons ( 



1991 



9841) . X ray flares (IRieger et al 



Dimitropoulou et al 



Droge et al. 



microwave flares ( 



1984 



Dennis 



1985 



Bai &: Sturrock 



7 ray 



19871 : 



20081 ). flares produ cing energetic inter planetary 



1990h. tvpe I 



Bogart fc Bai 



absent during solar cycle 22 ( 



1985 



and IV radio bursts ( Verma et al. 



Kile fc Clivei 



Kile fc Clivei 



1991 



Bai 



19911), and 



9911). However, this period icity was 



1992a 



Ozgiic fc AtaG,1994h . 



The periodicity has also been detected in indicators of solar activity (sunspot blocking 
function, sunspot areas, "active" sunspot groups, group sunspot numbers) which suggest 
that it is a ssociated preferentially with photospheric regions of compact magnetic flek 



struc t ures ( 



1990 



Lean fc Bruecknei 



Bouwei 



Ballester et al. 



1992; 



1989 



Lean 



Carbonell fc Ballestei 



1999 



Krivova fc Solanki 



99G: 



1992 



*ap et al. 



1990: 



Verma et al. 



Carbonell fc Ba. 



1992 



Oliver et al. 



lester 



1998; 



20021 ). Probably, the most important, and 



enigmatic, feature of the periodicity is that it appears during epochs of maximum activity 
and that it occurs in episodes of 1 to 3 years. 



Rabin et al.l (Il99ll ) performed a study of the magnetic flux variations during solar cycle 
21 which reveals the existence of quasi-periodic pulses or episodes of enhanced magnetic 
activity. The duration of the pulses is ~ 5 rotations during the years around maximum 
activity, the epoch in which the flare periodicity appears, and the comparison with magnetic 
fleld maps indicates that those pulses of activity co rrespond to the occurrence of complex 



active regions containing large sunspots ((Bai 



1987ah 



Ballester et al 



(12002 



20041 ) analyzed several data sets of, or strongly related to, 
photospheric magnetic flux to point out that the appearance of the near 160-day periodicity 
in different manifestations of solar activity during solar cycle 21 has its underlying cause in 



the appearance of the periodicity in the magnetic flux linked to regions of strong magnetic 
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field. They also showed that during solar cycle 22 the periodicity does not appear in the 
photospheric magnetic flux records and, as a consequence, the periodicity did not appear in 
other solar activity indicators, while during solar cycle 23 it appeared in the photospheric 
magnetic flux but not in other solar activity indicators. 

Several mech a nisms have been put forward in order to explain the existence of this 



periodicity. 



Wolfj (Il983l ) linked it to the interaction of 



rota t ing fea tures (active longitude 



bands) resulting from g- modes with I = 2 and / = 3. 



Bail (Il987bl ) suggested that the 



cause of this perio dicity must be a mech anism that causes active regions to be more flare 



productive. Later, 



Bai &: Sturrockl (119871 ) concluded that it cannot be du e to the int e ractio n 



of "hot spots", i.e. regions where flare activity is higher than elsewhere ( iBai 



1987a 



19881), 



rotat ing at different rates and that the cause must be a mechanism involving the whole 



Sun. 



Ichimoto et al. 



(1l985l ) suggested that it is related to the timesc a. 



e for storage and/or 



Bai fc Cliverl fll990[ ). taking into 



escape of magnetic flelds in the solar convection zone, 
account the possible intermittency of the periodicity, suggested that this behavior could be 
simulated with a damped, periodically forced non-linear oscillator, which shows periodic 



behav ior for some values of the parameters and chaotic behavior for other values. 



Wolfi 



(119921 ) argued that such periodicity can be understood in terms of the normal modes of 
oscillation of a nearly spherical, slowly rotating star, when two r-modes (inertial modes) 
couple with an interior g-mode beat. This suggestion seems to agree q ualitatively with 



the fa ct th at the periodicity is st ronger around the activity maximum. 



( ]199lh and 



Sturrock fc Bai 



Bai fc Sturrock 



(119921 ) proposed that the Sun contains a "clock" , modeled by an 



oblique rotator or oscillator, with a period of 25.8 days and suggested t hat the periodic ity 



of 154 days is just a subharmonic of that fundamental period. Later, 



Bai fc Sturrock 



( I1993I ) modified the earlier period to the value 25.50 days, but that model seems to be 



very c onstrained by helioseismological data about the rotation of the Sun's interior. 



Lou 



( I2OOOI ) suggested that such periodicities can be related to large-scale equatorially trapped 
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Rossby-type waves showing that, for typical solar parameters, the periods of these waves 
(with n = 1 and m even) are in good agreement with the observed ones. Moreover 



Lou 



( I2OOOI ) has also pointed out that such waves can give rise to detec t able f eatures, such 



Kuhn et al. 



(120001 ) have reported 



as surface elevations in the photosphere. Coincidently, 
observations made with MDI onboard SOHO and claim to have detected a regular structure 
of 100-m "hills", uniformly spaced over the surface of the Sun with a characteristic 



separation of 90,000 km. They suggest that t his structure is the su r 



Rossby waves, or r-modes oscillations. Finally, 



Dimitropoulou et al. 



ace manifestation of 



fl2008f ) have linked the 



found periodicities in d i fferen t classes (B, C, M, X) of X-ray flares with the theoretical 



Lou 



(l2000l ). pointing out that odd m periodicities are also frequent and 



periods derived by 
significant. 

On the other hand, most of the proposed mechanisms to explain solar flares, special 



the most energetic ones, accept as a prerequisite the emergence of magnetic flux (jPriest 



1990 



Forbes 



y 



I991I ) which, by reconnection with the ambient field, triggers the destabilization 



of active regions. Based on this mechanism, Carbonell & Ballester (1990, 1992) suggested 
that the periodic increase in the occurrence rate of energetic flares is related to a periodic 
emergence of magnetic flux through the photosphere. Later, Oliver et al. (1998) showed 
that during solar cycle 21 there was a perfect time correlation between the intervals o " 



occurr ence of the periodicity in sunspot areas and energetic flares, and 



Ballester et al 



(I2OO2I ) clearly pointed out that in cycle 21, and during the time interval in which the 
periodicity appeared, there was a perfect time and frequency coincidence between the 
impulses of high-energy flares and those corresponding to strong photospheric magnetic 
flux. The effi ciency of the reconnec tion mechanism depends on the geometry of the two 



flux systems (jGalsgaa: 



Zuccarello et al 
production. 



rd et al 



20071 ) and recent high resolution observations performed by 



20081 ) have confirmed the suitability of the mentioned mechanism for fiare 
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Emerged magnetic flux is probably connected to deeper regions, namely to the 
tachocline, which is a thin, transition layer between differentially rotating convection zone 
and rigidly rotating radiative envelope. The tachocline may preven t the spreading of the 



solar angular mornentum 



Gough et al. 



1998 



Gough 



rom 



: he convection zone to the interior (jSpiegel fc Zahn 



2007 



1992 



Garaud 



20071 ) and probably it is the place, where the 



large-scale magnetic field which governs the solar activity is generated/amplified. 

The observed periodicity of 155-160 days in the emerging flux is in the range of 
Rossby wave spectrum. Therefore, we suggest that the periodicity is connected to the 
Rossby wave activ i ty in the tacho cline. Rossby waves are well studied in the geophysical 

1987D. however, the prese nce of magnetic fields significantly 



context ( Gill 



1982 



Pedloskv 



modifies their dynamics ( IZaqarashvili et al 



2007 



20091 ). On the other hand, the differential 



rotation, which is inevitably present in the tachocline, may lead to the instability of 
particular harmonics of magnetic Rossby waves. It has been shown that the joint action 
of toroidal n iagnetic field and the differentia l rotation generally leads to tachocline 



instabilities ( 



2007 



Gilman fc Fox 



Gilman et al. 



1997 



Callv 



2003 



Dikpati fc Gilman 



2005 



Gilman fc Cally 



20071 ). However, the stability analysis usually has been performed in 
an inertial frame, which complicates to extract the information about unstable Rossby 
modes. Therefore, it is of paramount importance to perform the stability analysis in a 
rotating frame. Another important point is that the cons ideration of a rota t ing fr ame may 



tighten the stability criteria as it has been suggest ed by 



difference between the present analysis and that by 



Hughes fc Tobias (1200 ih . The 



Hughes fc Tobias! (120011 ) is the inclusion 



of rotation which allows us to obtain Rossby wave solutions. 

In this paper, we use a rotating spherical coordinate system to study the linear stability 
of magnetic Rossby waves in the solar tachocline taking into account the latitudinal 
differential rotation and the toroidal magnetic field. We perform a two dimensional 
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analysis, which can be follow ed in the future by more sophisticated shallow water 

Gilmaru 120001) . We first derive th e analy tical conditions of instability similar 



co nsiderations ( 
to iDahlburg et al. 



fll998h and 



Hughes fc Tobias! (1200 ll ). Then, we perform a d etailed 



stability analysis using Legendre polynomial expansions (ILonguet-Higgins 



19681 ) to obtain 



the spectrum of unstable harmonics of magnetic Rossby waves. 



2. Magnetic Rossby wave equations in the presence of differential rotation 

and the toroidal magnetic field 

Since the Rossby wave spectrum is clearly seen in the rotating frame, in the following 
we use a spherical coordinate system (r, 6', (p) rotating with the solar equator, where r is the 
radial coordinate, 9 is the co-latitude and is the longitude. 

The solar differential rotation law in general is 

n = Qo + ^i{0), (1) 

with 

ni{e) = -no{s2cos^e + s^cos'^e), (2) 

where f^o is the equatorial angular velocity, and S2, S4 are constant parameters determined 
by observations. 

Rossby waves are mainly polarized in the plane per pendicula r to gravity, then a 



19821 ). The two-dimensional 



two-dimensional {9, 0) analysis is a good approximation (iGilll 
analysis is also justified by Squire's theorem which states that for each unstable 3- 
dimensional disturba nce there is a corresponding unstable 2-dimensional disturbance with 



stronger growth rate (ISquire 



I933I). 



The magnetic field is predominantly toroidal, B = Se^, in the solar tachocline, and 
we take S = 5,^(6') sin 6', where is in general a function of co-latitude. Then, the 
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incompressible magnetohydrodynamic (MHD) equations in the frame rotating with Qq are 
(see appendix A): 



due 
dt 



dufj 



2[no + ni{9)] cos 9u^ 



1 dpt 



■ — -i— ;^ — Or/ 



4:7ipRo 



^ + fil (0) ^ + 2^0 cos Que + fii (^) cos Oue + ue-^ [sin ^fii (^)] 
at od) oU 



Rq sin 9 



dpt _^ B^ db^ _^ 



be 



d 



A-jipRo 



ATrpRo sin 9 89 



{B^ sin' 



B^due d db^ 



-(smM + ^ = 0, 



(3) 



(4) 
(5) 
(6) 



where ug, u^p, bg and b^ are the velocity and magnetic field perturbations, pt is the total 
pressure (hydrodynamic plus magnetic), p is the density and -Rq is the distance from the 
solar center to the tachocline. 



We consider the stream functions for velocity and magnetic field 



Ug 



1 9^ 

sin 9 dd) 



u^ 



5^ 



1 9$ 
sin 9 dd) 



9$ 



(7) 



Substitution of expressions ([7]) into ©-([SD and Fourier analysis with exp[im(0 — ct)] 



gives 



d 



d 



09^^^^ 89 sin 9 



d 



^ - 2^]o sin + — 



1 d 



d9 V sin 9 d9 



{VLi siv? 



B, 



4:7ipRo 



d . „ d m? 



d 



1 d 



(c-fii)$ = -^vl/. 



4:7!- pRod9 \sm 9 d9 

B, 



{B^ sin^ 



(9) 



Let us now make the transformation of variables p = cos 9, then we obtain (\E' and $ 
are normalized by fio-Ro and Bq respectively, where Bo is the value of B^ at 9 = 0) 

(Qd - i^)L^ + (2 - - - I3^BL<^ + P^-^[B{1 - P^)]<^ = (10) 
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where 

is the Legendre operator and 



d_ 



d 



m 



9/i 1 — /^^ 



UJ 



Bl 



Eqs. ( ITO|) - f|TT|) govern the 2- dimensional dynamics of magnetic Rossby waves in the 
presence of differ ential rotation and tor oidal magnetic field. The equations are analogous to 



Eqs. (17)-(18) of iGilman fc Foxl (119971 ). but are written in the rotating frame instead of in 
the inertial one. 



3. Analytical conditions of magnetic Rossby wave instability 



In tl 


lis section, we derive t 


le analytical instability bounds using a well known technique 


( Howard 


1961; 


Drazin & Reid 


1981; 


Watson 


1981; 


Gilman &; Fox 


1997; 


Dahlbure et al. 




1998; 


Hughes & Tobias 


2001 


)• 





Let us define a new function H 

^ = - u)H, $ = BH. 
Then Eqs. f[TO|) - f|TT|) can be cast in the following form 

|-(1 - - j^P{fi)H + 2{n, -u;)[l + {f,n,y]H - 2p'B{pByH = 0, (12) 

where 

P{^i) = {Q,-uj)'-p'B' 
and ' means differentiation with respect to fi. 
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Now, multiplying Eq. f|T2|) by H*, integrating from -1 to 1 and using the boundary 
conditions H{ii = ±1) = 0, we get 



j P{fi)Qdfi- j 2{nd-uj)[l + {nndy]\H\'^dn + j 2P^B{iiBy\H\'^dfi = 0, (13) 



where 



Q = {l-fi' 



dH 



+ 



m 



1-/.2 



\H\' > 0. 



Considering u = Ur + ioOi in Eq. (fT3|) we obtain two different conditions for instability 
(see detailed derivations in appendix B). The first condition states that the instability takes 
place when 



(14) 



with 

fif=[(.,f.^ + .4/)^-,3V]_. 
In the remaining ^ax and mean maximal and minimal values. 



(15) 



This means that the frequencies of unstable harmonics (actually phase speeds, while 
frequencies can be obtained by multiplying by m) lay inside the upper semicircle of complex 
cij-plane with center at the origin and radius Ri (see Fig. 1). 



1 r 



The second instability condition is the semicircle theorem similar to iHoward fil96ll ). 



The MHD ge neralization of Howard 's se micircle theorem in rectangular coordinates has 



been done by 



Dahlburg et al 



mm and 



Hughes fc Tobias! (l200ll ). Here the theorem is 



derived in the rotating spherical coordinate system as the second condition of instability 
(see details in appendix B), obtaining 
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2 / \ 2 

^dmin H" ^dmax \ , 2 / ^dmin H" ^dmax 



where 



1 -/i^ 



'] + ^-^2P^B{fxBy - P^B^. (17) 



We observe that fl 



dmax 



and 



dmin 



-e, where e = S2 + S4, therefore we can write 



+ -j <— + Amax- 



Due to this condition the frequencies of unstable modes lay inside the semicircle of the 
complex cij-plane with center 

(19) 



and radius (see Figure 1) 



R-, 



(20) 



Equations flMj) and f[T8|l are two necessary conditions of instability. They define two 
different semicircles in the complex a;-plane, an d the instability occurs when the two 



semicircles overlap (see 



Hughes fc Tobias! (120011 ) for the same statement in the rectangular 



case). If the radius of one semicircle tends to zero, the instability disappears. 



In the remaining we use a magnetic field 



B4, = Bqh, 



(21) 



which changes sign at the equator (IGilman fc Foxl 119971 ). 



Now, we may estimate the instability bounds under tachocline conditions. An 
important step is to choose the parameters of differential rotation, S2 and S4. These 
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parameters are determined by observations and their values at the solar surface are 
S2 ~ S4 ~ 0.14. Helioseismology shows that the transition between the differentially 
rotating convective zone and the rigidly rotating radiative interior is described by the 
function $(r, Tc, w) = 0.5(1 + er/[2(r — rc)/w]), where erf is the error function, Vc is the 
radius of the central point of the tachocline and w is the characteristic thickness of the 
tachocline corresponding to a variation of ^(r) from 0.08, at the bottom of the tachocline. 



to 0.92, at the tachocline's upper surface (iKosovichev 



19961 ). In order to calculate the 



parameters of the differential rotation at the upper part of the tachocline, the solar surface 
values must be multiplied by 0.92, then, we obtain S2 ~ S4 ~ 0.13. However, it must 
be mentioned, that the re al values of these parameters can be different in the tachocline 



( ICharbonneau et al. 



oscil 



2000 



ations ( 



Howe 



19991) and also can change through the solar cycle due to torsional 



LaBonte fc Howard 



1982 



Komm et al. 



1993 



Antia fc Basu 



200c 



Howe et al. 



2OO9I ). Therefore, these values are tentative and further observations are needed 



to infer the correct parameters and their cycle dependence. 

The typical values of equatorial angular velocity, radius and density in the tachocline 
are Qq = 2.7 ■ 10^^ s^^, Rq = 5 ■ 10^° cm and p = 0.2 g ■ cm~'^ respectively. Then, the 
parameter is much smaller than unity being ^ 0.0022 for a magnetic field strength of 10^ 
G. Using these parameters we get -Ri = 0.256 and R2 = 0.154 for azimuthal wave number 
171 = 1. Then, the conditions f|T^ and f|T8l) give that the minimum period of the m = 1 
unstable modes in the tachocline is 



T 



105 days. 



(22) 



Therefore, only the magnetic Rossby modes with periods longer than 105 days may grow 
in time. However, equation fl22]) only gives a lower bound for oscillation periods. A more 
detailed analysis is required to reveal the spectrum of possible unstable harmonics. 
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4. Spectrum of unstable magnetic Rossby modes 

In this section, we use the general technique of Legendre polynomial expansion 



( iLonguet-Higginsi Il968l ) . Using the magnetic field profile (|2T1) . Eqs. f|TOl) - f|TTl) are rewritten 



as 



(rid - + (2 



rf2 



(23) 
(24) 



Let us expand ^ and $ in infinite series of associated Legendre polynomials 

cx;' oo 

VI/ = 5^ a„pr(/x), $ = 5^ &nPr(/i), 

n=m n—m 

which satisfy the boundary conditions \E' = $ = Oat/i = ±1. 

The latitude-dependent part of the differential rotation has the form 



(25) 



(26) 



We substitute (!25|) into Eqs. (!23l) - (!2^ and, using a recurrence relation of Legendre 
polynomials, we obtain algebraic equations as infinite series (Details of the calculations can 
be found in Appendix C for the case when the differential rotation has only second order 
dependence on fi in expression fl26|) ). The dispersion relation for the infinite number of 
harmonics can be obtained when the infinite determinant of the system is set to zero. In 
order to solve the determinant, we truncate the series at n = 75 and solve the resulting 
polynomial in u numerically. The frequencies of different harmonics can be real or complex 
giving the stable or unstable character of a particular harmonic. It turns out that m = 1 



in many differen 


t occasions ( 


Watson 


1981; 


Oilman & Fox 


1997; 


DikDati & Oilman 


2005 


Oilman & Callv 


2007 


)■ 
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Figure [2] shows the real, mcr, and imaginary, mcj, frequencies of all m = 1 unstable 
harmonics for different combinations of differential rotation parameters and magnetic 
field strength. In order to show the dependence on the parameters ^2,54, we vary these 
parameters for different values of magnetic field strength so that the sum S2 + S4 (which 
is the difference in equatorial and polar angular velocities) remains 0.26. In Figure O 
the upper left panel corresponds to the case considered in Appendix C (i.e. S4 = 0). 
Blue, green, yellow and red colors correspond to magnetic field strengths of 2 ■ 10^ G, 
6 ■ 10'^ G, 2 ■ 10^ G and 4 ■ 10^ G, respectively. Asterisks (circles) denote the symmetric 
(antisymmetric) harmonics with respect to the equator. The results show that the 54/1^ term 
i n the differential rotation (l26l) significantly affects the behaviour of unstable harmonics 



(ICharbonneau et al. 



I999I ). For each combination of S2, S4 and the magnetic field strength, 
there is a particular unstable harmonic with a growth rate much stronger than for the other 
harmonics. This harmonic is symmetric with respect to the equator and has the frequency 
of 0.17-0.18 Qq (yielding periods of 150-160 days) for the magnetic field strength of < 2 ■ 10^ 
G. The frequency decreases for stronger magnetic fields (red colors), therefore Rieger-type 
periodicities arise as symmetric unstable harmonics for relatively weaker magnetic field 
strength. 

Thus, the appearance of a strong oscillation with a particular frequency needs a 
suitable combination of differential rotation parameters {s2, S4) and magnetic field strength. 
However, the differential rotation parameters used in Figure |2] are probably too high for 
the solar tachocline. Therefore, we study the dependence of unstable harmonics on more 
realistic differential rotation rates. 

Figure [3] displays the dependence of the most unstable symmetric harmonic (this 
harmonic can be identified on Fig. 2 as the blue, green, yellow and red asterisks at top of 
each panel) on the differential rotation parameters for two different values of the magnetic 
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field. Left panels correspond to tlie field strength of 2-10^ G and right panels correspond 
to the strength of 10'^ G. Real and imaginary parts of the harmonic vs S4 are plotted for 
different values of S2- The values of S2 vary from 0.14 (blue color) to 0.09 (yellow color). 
We can observe that the frequency, mcr, of this harmonic is only slightly dependent on the 
differential rotation parameters and takes values between 0.16-0.18 Qq which correspond to 
oscillation periods of 150-170 days. This is the range where the Rieger-type periodicity has 
been observed. On the contrary, the growth rate, mcj, of this harmonic strongly depends 
on the differential rotation parameters. The growth rate becomes stronger when both S2 
and S4, are increased. 

The frequency and growth rate of this harmonic have no significant dependence on the 
magnetic field when its strength is smaller than 10^ G. Figure H] shows the dependence 
of the harmonic calculated for three different profiles of the differential rotation (blue line 
corresponds to S2 = 0.13, S4 = 0.1; the red line to S2 = 0.11, S4 = 0.12 and green line to 
S2 = 0.11, S4 = 0.1). We can observe that the stronger growth rate occurs for the red line, 
which means that S4 is more important for the instability. 

When the magnetic energy becomes comparable to the energy of differential rotation, 
then the frequency of the symmetric harmonic is significantly reduced (see red asterisks on 
Figure [2]). The critical magnetic field strength, i.e. when the magnetic energy is comparable 
to the flow energy, is ~5 ■ 10^ G for the differential rotation parameters S2, S4 =0.13. In 
this case, (s2 + ^4)^ ~ the radius of first semicircle Ri (see Eq. (15)) tends to zero and 
the growth of symmetric unstable harmonics is suppressed. 
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5. Discussion 

The periodicity of 155-160 days was discovered almost three decades ago, however 
the reason of its appearance/disappearance is still unknown. The most striking feature, 
perhaps, is its appearance only at certain times, which normally coincide with the maximum 
of the cycle (Figure [5]). This coincidence naturally suggests that the magnetic field and 
the differential rotation at the solar cycle maximum provide suitable conditions for the 
appearance of this periodicity. 

Here we show that the periodicity can be connected to the dynamics of magnetic Rossby 
waves in the tachocline, since, in this layer, they are unstable due to the presence of toroidal 
magnetic field and latitudinal differential rotation. First, we have derived the analytical 
bounds of instabihty, which state that m = 1 unstable modes have periods > 105 days. 
Next, we have calculated the detailed spectrum of unstable harmonics using the method of 
Legendre polynomial expansion. We have found that the behaviour of unstable harmonics 
is very sensitive to the combination of magnetic field strength and the differential rotation 
parameters (52,54). Each combination of the parameters favours a particular harmonic, 
which has stronger growth rate compared to other unstable harmonics. Therefore, this 
harmonic may quickly dominate over the others and may lead to a detectable oscillation, if 
the parameters remain more or less unchanged during some time. Unstable harmonics have 
two types of symmetry with respect to the equator: symmetric and antisymmetric. The 
growth rates of symmetric modes are higher than the antisymmetric ones and they depend 
on the differential rotation parameters; the growth becomes stronger for stronger shear. 

Frequencies of symmetric unstable modes are in the range 0.16-0.18 VLq (Figure [3]), 
which yield the periods of 150-170 days. In the case of strong differential rotation, their 
growth rate may reach up to 0.015 VLq i.e. the growth time is ~ 280 days. Therefore, they 
may quickly dominate over the rest. The growth of the magnetic Rossby wave amplitude 
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leads to an enhanced magnetic buoyancy at the tachochne which causes the periodic 
eruption of magnetic flux towards the solar surface. Therefore, the periodicity is observed 
in the emerged magnetic flux and consequently in many indicators of solar activity (see 
references in the Introduction). 

The question why the periodicities appear only at particular times (mostly just after 
solar maximum, see Figure [5]) needs additional explanation. A possible reason is that the 
growth of symmetric harmonics strongly depends on the differential rotation parameters 
(32,34). It is known that the solar differential rotation is changing through the solar 
cycle. The pattern known as the torsional oscillation has been first observed at the solar 
surface in full disc velocit y measurements (U 



magnetic features as well ( iKomm et al 



LaBonte fc Howard 



19821 ) and later in surface 



I993I ). Helioseismology shows that the torsional 



oscillatio n is not only a sur: 



interior (lAntia &: Basu 



2000 



ace phenomenon 



Howe et al 



3 ut may pen etrate deeper into the solar 



2000 



Howe 



20091 ) . Then, the parameters ^2,^4 



may vary through the solar cycle in the tachocline, which permits the strong growth of 
symmetric magnetic Rossby waves only at particular times. This time should coincide with 
the solar maximum. We think that additional helioseismic estimations are needed to study 
this phenomenon. 

One of the significant simplifications in our approach is the linear stability analysis. 
The growth of perturbation amplitudes probably leads to nonlinear effects. On the other 
hand, the process would be accompanied by increased magnetic buoyancy, which causes the 
eruption of magnetic flux upwards and consequently may stop further growth of amplitudes. 
These processes should be studied with sophisticated numerical simulations in the future. 



It shoulc 
instabilities ( 


be mentioned 


lere t 


hat numerous previous papers I 


lave studied the tachoc 


ine 


Gilma 


n & Fox 


1997; 


Gaily 


2003 




Dikuati & Oilman 


2005 




Oilman & Oallv 


2007; 


Oilman et al. 


2007 


). However, all the calculations have been performed in an inertial 
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frame, while the Rossby wave dynamics is more clearly seen in a rotating frame. Another 
important difference between inertial and rotating frar nes is that the instability conditions 
may be tightened in the moving frame as suggested by lHughes &: Tobias! (120011 ) . 



The solar tachocline may consists of two parts: the inner radiative layer with a strongly 



stabl e stratification and the outer overshoot layer with a weakly stable stratification (IGilman 
200ol ). The latitudinal differential rotation should be stronger in the upper tachocline and 
weaker in the lower one. On the contrary, the magnetic field strength should be higher in 
the lower part and smaller in the upper one. Therefore, the upper tachocline may favor 
the better conditions for the growth of symmetric unstable harmonics, which trigger the 
Rieger-type periodicities. 



6. Conclusions 

In summary, we have shown that the destabilization of magnetic Rossby waves in the 
solar tachocline is produced by the joint effect of the latitudinal differential rotation and 
the toroidal magnetic field. The frequencies and growth rates of unstable harmonics depend 
on the combination of the differential rotation parameters and the magnetic field strength. 
The possible increase of latitudinal differential rotation at the solar maximum may trigger 
the instability of symmetric harmonic with period of 155-160 days in the upper part of the 
tachocline. This instability has a direct correlation with magnetic flux emergence, therefore 
the periodicity also appears in solar activity indicators related with magnetic flux. Later on, 
and probably via reconnection, this periodic magnetic flux emergence triggers the observed 
periodicity in solar flares. 

The magnetic Rossby wave theory opens a new research area about the activity on 
the Sun and other stars, and magnetic Rossby waves can be of paramount importance for 
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observed intermediate periodicities in solar and stellar activity (IMassi et al.lll998l . l2005l ) 
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A. MHD equations in a rotating frame 

In the case of rigid rotation it is straightforward to transform equations from inertial 
into the rotational frame, but the presence of differential rotation slightly complicates the 
considerations as different parts of the Sun rotate with different angular velocity. The best 
way to overcome the difficulty is to consider the frame rotating with the equator. Then 
the latitudinal differential rotation can be considered as the unperturbed shearing motion 
in this frame. 2-dimensional incompressible linearised MHD Equations {6, 0-plane) in the 
frame rotating with angular velocity of the equator, Qq, are 

due , due oo a o cos6' 1 % , H dbe ^ H cos 61 

— — 2\l()COSauj,—2 U^u^ = — 2 Oa 

dt RosmO " * RosinO ^ pRo d9 4n pRo sin 9 d(f> ATrpRosmO^ 

(Al) 

du^ du^ ugdU^, cose 

-m + R^^e^ + + + R^e""'^' = 

' ^'''+ + (A2) 



Rq sin 9 (90 AnpRo sin 9 dcj) AnpRo sin 9 89 



^The software is available at http://paos.colorado.edu/research/wavelets 
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dbe ^ U4, dbg 
dt Rq sin 9 d(j) 



E due 



Rq sin 9 
|(sinM + ^-0, 



(A3) 
(A4) 
(A5) 



where ug, u^p, bg and b^ are the velocity and magnetic field perturbations, S and are 
azimuthal components of unperturbed magnetic field and velocity in the rotating frame, pt 
is the perturbation in total (hydrodynamic plus magnetic) pressure. 

We consider [/^ as the differential rotation with respect to the equator, i.e. 



= Rosm9ni{9). 



(A6) 



Substitution of this expression into Eqs. (lAip - (lA5p gives Eqs. 



B. Derivation of analytical instability conditions 



The real and imaginary parts of Eq. fll3l) with u = Ur + iuji are 



and 



2iuJi 



{nd-ujr)Qdfi- [1 + {findy]\H\'di2 



0. 



Unstable harmonics should have non zero Ui, therefore Eq. (IB2p requires 



-1 



:bi) 



:b21 



Substitution of J^^ ^IdQdfi from this equation into Eq. (IBIO leads to the equation 



- 21 - 



(B3) 



which then can be rewritten as 

"1 



1 



P'B'] (1 - 



OH 



+ 2p^B{iiB) 



m 



This equation will be satisfied if both integrals are zero, which requires 



\H\^dyi = 0. 



{nl - I3''B-'U^ < ujI + ujI < {nl - (B4) 

and 



Q2 _ ^2^2 _ 2^^[^ ^ (^^^)']L^ + 2f3'B{f^By^—^ ) < + uj^ < 



mm 



Ql - (3'B' - 2Q,[1 + (^Q,y]L^ + 2(3'B{^,By^—^ ) . (B5) 



max 



The inequality ( ]B5ll is similar to the inequality (1B4I) . but with two additional terms 
in the left and right hand sides. Both additional terms are positive, therefore, inequality 
( ]B4p determines a condition of instability. Using the profiles of magnetic field ( !2TI) and the 
differential rotation ([2]), Eq. (IB4p leads to Eq. (ITll) in the main text. 

In order to obtain the semicircle theorem let us observe that 

J {^Id - ^ldmin)i^d - ^dmax)Qdll < 0. (B6) 

Then the substitution of j^-^ Vt^^Qd^ from Eq. (1B3P into Eq. (1B6P gives 

j \ujI + + 0^B^ — {Vtdmin + ^dmaz)l^r + ^dmin^dmax\ QdfX < 
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This inequality can be rewritten as 



+ tUj + p i3 — + i Idmin}' ^dmax 





OH 







2 



J V ^ 2 y H-CUj +p -D — I I + iZdmm"dmaa; — 

1 — U'^ 1 — Zi'^ ?7i^ 

— ^(n^™,, + n^^ax - + (i^n,)'] 1^2/325 (/iB)']- — ^i/^i^d/i < 0. (B7) 

m"^ 1 — /i"^ 

At least, one of the two integrals should have negative sign, therefore 



""arrmi i ""arnax \ 2 , / n2 t:j2\ f ^dmin ~^ ^dmax\ ^ n /'DO^ 
+ l^i + IP -D )min — I + ^ l-dmin^ Umax S U {tSo} 



and / or 



^ I ~r I ^ I ~r ^ I'dmin^ ^dmax ^max _ U, l-t^yj 



where 



A{ii) = ^-^{n,^,^ + n,^ax - 2f^d)[l + {iin^)'] + ^-!^20'B{iiBy - 0'B\ (BIO) 



The inequality flB9P is wider than (IBSp . Therefore, it determines the second condition 
of instability (Eq. flTBl) in the main text). 

C. Derivation of dispersion equations using Legendre polynomial expansion 



The substitution of (|25|) into Eqs. fl23|) - fl24l) and using the Legendre equation 
LP^ + n(n + 1)P™ = leads to 



°° j2 
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- ^ bnPJT = 0, (CI) 

n=m 

oo oo 

{n, - ^) ^ 6„pr = E ^-^n- (C2) 

n=m n=m 

Now we take the explicit form of Qd = —^fJ''^, then 

oo oo oo 

J2[^n{n + l) + 2 + 2e]a„Pf + e ^ + 1) - 12]a„/i2pf + /J^ ^ + 1) - 6]6„/iP„"^ = 0, 

n=m n=m n=m 

(C3) 

oo oo oo 

J2 ^nliP'n + 5^ c.6„Pr + e 5^ 6„/i2pr = 0. (C4) 

n=m n—m n—m 

We use the recurrence relations between Legendre polynomials, namely: 

2pm _ A pm I n pm , ^ pm 
,. pm _ p) pm I p pm 



where 



^ (n + m)(n + m — 1) ^ {n — m){n + m) {n — m + l){n + m + 1) 
~ (2n + l)(2n-l) ' " ~ (2n + l)(2n- 1) ^ (2n+ l)(2n + 3) ^ 

(n — m+l)(n — m + 2) ^ n + m ^ n — m + 1 
"~ (2n+ l)(2n + 3) ' " ~ 2n + 1 ' " ~ 2n + 1 ' 

Substitution of these relations into Eqs. (1C3I) - (1C4I) gives 



[^n{n + 1) + 2 + 2e]a„Pr + e + 1) " 12] A„a„Pr_2 + e J] + 1) - 12]P„a„Pr + 

n—m n=m n=m 

oo oo oo 

+e 5^[n(n+l)-12]aa„P„";.2+/3' 5Zh(n+l)-6]D„6„Pr_i+/52 5^[n(n+l)-6]E„^^^ = 0, 

n=m n=m n=m 

oo oo oo oo oo oo 

n=m n=m n=m n=m n=m n—m 

Rearranging terms we obtain 
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5^[a;n(n+l)+2+2e]a„Pr+e$^[(«+2)(r^+3)-12]A„+2a„+2Pr+e5Zh(^+l)-12]5„a„P„^ 

n=m n=m n=m 

oo oo 

+e [(^ - 2) (n - 1) - 12]C„_2a„_2Pr + J] [(^ + l)(n + 2) - 6]6„+iZ}„+iPr+ 

n=m n=m 

oo 

J2 Hn - 1) - 6]6„_iE„„iP™ = 0, 

n=rra 

oo oo oo oo oo 

an+l/^n+lP^ + 5^ an^iEn-lPT + ^^nP^ + ^ A„+2&n+2Pr + ^ BnhnPn + 

n=m n=m n=m n=m n=m 

oo 

+e ^ C„_2&„-2Pr = 0. 

n=r?i 

Now the coefficients of P™ give the equations 

0, (C5) 
+ P„6„+2 + M„6„_2 = 0, (C6) 

where 

, , r / N — m)(n + m) , , , , (n — m + + m + 1) 

S„ = .n(n+l)+2+2.+.[n(n+l)-12] ^'^^^ ^ ^^^^^^ _ +.ln(n+l)-12| ^ ^ 

An + m + 2)(n + m + 1) ^ ,, , Jn — m — l)(n — m) 

F„ = .|(„+2)(„+3)-121 ^ (2„ + 5)(2„ + 3) ' ^» = ^l'-^"-^'-"! (2n - 3)(2; - 1) - 

n — m [n — m){n + m) {n — m + l){n + m + 1) 

""^2^^' ^""'^ + '(2n + l)(2n-l) (2n+l)(2n + 3) ' 

^ (n + m + 2)(r;, + m + 1) ^ (n — m — l)(n — m) 
(2n + 5)(2n + 3) ' " ~ ^ (2n-3)(2ra-l) 

The expressions ( ]C5l) -( !C6l) are infinite series and the dispersion relation for the infinite 
number of harmonics can be obtained when the infinite determinant of the system is zero. 
In order to solve the determinant, we cut the series at n = 75 and solve the resulting 
polynomial in u numerically. 
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Fig. 1. — Semicircles of unstable harmonics in the complex {cUr, a;i)-plane corresponding 
to the two instability conditions, Eqs. (14) and (18). Instability occurs when these two 
semicircles overlap, ujr, Ri and R2 are normalised with respect to f^o- 
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Fig. 2. — Real [mcr) vs imaginary (mcj) parts of unstable harmonic frequencies for different 
combinations of differential rotation parameters S2, S4 and magnetic field strengths (frequency 
is normalized by equatorial angular velocity, Qq)- Note that the difference between equatorial 
and polar angular velocities S2 + S4 = 0.26 remains the same for all panels. The toroidal wave 
number m equals 1. Blue, green, yellow and red colors correspond to magnetic field strengths 
of 2 ■ 10^ G, 6 ■ 10'^ G, 2 • 10^ G and 4 ■ 10^ G, respectively. Asterisks denote the symmetric 
harmonics with respect to the equator, while circles denote the antisymmetric ones. The 
frequencies are normalized by equatorial angular velocity, Qq, for example, mCr = 0.18 
corresponds to the period of ~ 150 days. 
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Fig. 3. — Real (lower panels) and imaginary (upper panels) part of the frequency of the 
most unstable symmetric harmonic vs S4 for different values of S2- Dark blue, green, red, 
blue, magenta and yellow colors correspond to 0.14, 0.13, 0.12, 0.11, 010 and 0.09 S2 values 
respectively. The magnetic field strength equals to 2-10^ G (left panels) and 10^ G (right 
panel) respectively. 




Fig. 4. — Dependence of real (lower panel) and imaginary (upper panel) part of the frequency 
of the most unstable symmetric harmonic on the magnetic field strength for three different 
combination of differential rotation parameters. The blue, green and red lines correspond to 
(s2=0. 13,54=0.1), (s2=0.11, S4=0.1) and (s2=0.11, S4=0.12) respectively. 
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Fig. 5. — Top Panel: Plot of the daily (black) and monthly averaged (red) sunspot areas 
for solar cycles 19 to 23. Bottom Panel: Time/period diagram calculated using the Morlet 



wavelet (ITorrence fc Compo 



19981 ) with fco = 20. Vertical solid white lines mark the epochs 
of minimum solar activity, while the two dashed lines correspond to the maximum of cycles 
19 and 21. Large power values around 160 days can be seen in cycles 19, 20, 21 and 23, 
peaking at cycle 19. Power is given in arbitrary units. 



